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Abstract 

We study percolation on the worldsheets of string theory for c = 0, 1/2, 1 and 
2. For c < 1 we find that critical exponents measured from simulations agree 
quite well with the theoretical values. For c = 1 we show how log corrections 
determined from the exact solution reconcile numerical results with the KPZ 
predictions. We extend this analysis to the large c regime and estimate how 
finite-size effects will effectively raise the ground state energy, masking the pres- 
ence of the tachyon for moderate values of c > 1. It thus appears likely that 
simulations for c = 2,3... on numerically accessible lattices will fail to even 
capture the qualitative behavior of the continuum limit. 
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1 Introduction 



Within the last few years, the critical behavior of low-dimensional (c < 1) string 
theories has been understood in considerable detail. In particular, both analytic 
and numerical work has yielded a description of the intrinsic geometry of the 
worldsheet of these models [1, 2, 3, 4]. The worldsheet has been shown to be 
very spiky and branched; it can be characterized by a scaling distribution of 
embryonic universes [5, 6]. As c increases towards 1, embryonic universes of all 
sizes but with necks of the order of the cutoff proliferate and perhaps saturate 
the worldsheet. 

A direct determination of properties of the intrinsic geometry of string the- 
ories constructed from c > 1 matter has proven much more difficult. There are 
many theoretical indications that the worldsheet degenerates into polymer-like 
configurations above c = 1. Though the c < 1 surfaces are quite spiky, they still 
appear to retain some of the character of two dimensional surfaces; this perhaps 
is no longer true for c > 1. For at c = 1 the dressed identity operator (coupled 
to the cosmological constant) becomes tachyonic. Above c = 1, this operator 
should then spawn states which are non-local in the worldsheet metric that effec- 
tively tear the worldsheet apart [7, 8] . Related computations in Liouville theory 
indicate that a sort of Kosterlitz Thouless transition is anticipated at c = 1, in 
which vortex configurations of the Liouville field become unbound [9] . Unfor- 
tunately though, this picture has not been verified via an exact calculation for 
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c > 1; the appropriate matrix models have not been amenable to exact solution 
(though see [10] and [11]) and Liouville theory yields complex exponents in this 
regime. We would therefore like to detect, through Monte Carlo simulations, 
some characteristic of the geometry of worldsheets that distinguishes c < 1 from 
c > 1 and signals the onset of the tachyon. 

In fact, numerical work has not provided any evidence of a dramatic change 
in the internal geometry of the worldsheet as c passes 1. Simulations, albeit on 
modest-sized lattices, only have shown convincing evidence of branched-polymer 
structure on the worldsheet for c > 10 — 12. Recently, there has been much 
work on simulations of multiple Potts models coupled to gravity [12, 13, 14]; 
the systems studied generally have small values of c both greater and less than 
1. In these studies, the internal geometry was characterized by measuring the 
distribution of discrete curvature, or more precisely, the distribution of ring 
lengths on </> 3 graphs dual to triangulations. These measurements indicated 
that no dramatic change in these curvature distributions occurs as c is increased 
above 1. Interpreting these results rigorously is difficult, though, sinc^rj is not 
clear how to identify these curvature distributions with continuum ( uni versal) 
correlation functions. 

To search for a transition to branched polymers, Ambj0rn et al.lZph] also 
simulated Ising models on typical fixed random triangulations sampled from the 
set of worldsheets characteristic of gravity coupled to matter of central charge 
c. Each triangulation was generated during a simulation of a multiple Potts or 
Gaussian models coupled to gravity; the Ising model was then simulated on each 
of these particular triangulations. These authors showed that the Ising model on 
a branched polymer will not undergo a phase transition. Only when c reached or 
exceeded 12 did they find that no finite-temperature Ising transition was evident 
on the surfaces extracted from their simulations. 

For fairly low central charge, the coupling between matter and gravity, as 
inferred from the above numerical work, appears to be fairly weak, perhaps too 
weak to drive the worldsheet into a branched polymer phase [12]. It has thus 
been suggested that worldsheets of gravity coupled to matter with c > 1 but less 
than about 10 might lie in some sort of intermediate phase and not undergo a 
transition to branched polymers for c above 10 or so [15, 12]. Such a scenario is 
quite intriguing, especially since the existence of such a phase does not seem to 
be anticipated by Liouville theory. 

To reveal the character of the critical geometry of string worldsheets, we 
shall examine site percolation on dynamical triangulations coupled to matter. 
The distribution of percolation clusters shouldrSeTve as a sensitive probe of the 
branched structure of the worldsheet, distinguishing c < 1 from a branched poly- 
mer phase. Essentially, embryonic universes act as traps that prevent percolation 
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clusters from growing, while regions of high connectivity will enhance the span 
of percolating clusters. We also indicate how the critical behavior of ^erjcolating 
clusters on worldsheets of c < 1 matter coupled to gravity can be determined 
from Liouville theory. The quality of the results of numerical simulations of 
two-dimensional gravity has often been called into question, with the suspicion 
that the lattices used were often too small too even qualitatively reflect contin- 
uum behavior. We shall see however that we can reproduce on lattices of several 
thousand nodes, with quite good precision (of order 1%), the theoretical predic- 
tions of these percolation exponents for c < 1. At c = 1, matrix model solutions 
and Liouville theory [16] predict a more complicated logarithmic dependence of 
correlation functions on the cosmological constant. We also find that our results 
are consistent with the functional dependence predicted by matrix models along 
with the exponents of Liouville theory. Our data shows only an apparent weak 
coupling of c = 1 matter to gravity, but this is misleading. The logarithmic 
scaling of the matrix model implies that effectively c = 1 matter only has a 
small effect on the percolative structure of the lattice for small lattice sizes. For 
larger lattices however, the matrix model solutions predict that the influence of 
the matt^F^oupling to gravity should gradually become more pronounced. 

We also will demonstrate that percolation on non-interacting branched poly- 
mers does not undergo a transition for p < 1. We shall then find evidence sug- 
gesting that for c = 2, the critical behaxklr of percolating clusters (or indeed 
perhaps the lack of criticality) may be qualitatively different from that of c < 1. 
The c = 2 observations, which need to be interpreted with great care, at first 
glance appear to be indicative of a breakdown of behavior that characterizes 
percolation on surfaces. 

In illuminating the above results, we shall argue that for c somewhat greater 
than one, the finite size of the lattice will largely mask the theoretically an- 
ticipated degeneration of the worldsheet. The infrared cutoff should induce a 
large shift in the tachyonic energy, preventing the unrestrained proliferation of 
tachyons on not too large lattices. We can, in fact, estimate the magnitude 
of this effect. With this in mind, we then argue that for c somewhat greater 
than 1, we are most likely measuring finite-size artifacts; yet the presence and 
qualitative nature of these effects can be understood in the context of string 
theory. 

A delicate analysis of percolation on dynamically triangulated lattices is 
possible largely due to a curious but elementary graphical property of trian- 
gulations: they are self-matching. This constrains the value of the percolation 
threshold to be 1/2, with a few caveats. We shall explain this further in Section 
3. 

In a companion paper [17], an analysis of the critical properties of the 
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two-species Ising model coupled to two-dimensional gravity appears. As in this 
paper, it is essential to incorporate logarithr riic b orrcctions to scaling to compare 
numerical results with theory. Percolation is also simulated on random lattices 
coupled to two Ising species, yielding results consistent with those preseniicl 
here. 

The plan of this paper is as follows. We begin by reviewing basic facts 
about percolation and discuss how the KPZ [18] equation predicts the behavior of 
percolating clusters on c < 1 dynamical lattices. We then explain the theoretical 
predictions for c = 1 and also analyze the behavior of percolation on branched 
polymers. A discussion of the constraint p c = 1/2 follows. We then describe the 
numerical techniques we have used, including an analysis of auto-correlations. 
We next present the results of our simulations (for c = 0, 1/2, 1 and 2) and then 
give our conclusions. 

2 Theoretical Predictions 

We shall present the results of site percolation on random lattices sampled from 
simulations of quantum gravity coupled to matter. Site percolation describes the 
following process: we randomly color sites 'black' on a lattice with probability 
p. Adjacent black sites are then connected together to form clusters. We shall 
generally be interested in features of the distribution n(s) of these clusters as 
a function of the number of constituent sites s. Typically, we will measure the 
mean cluster size, denoted by S. If one randomly chooses a black point on the 
lattice, then on average it will belong to a cluster of size S. Explicitly, S = 
(s 2 )/(s), in which averages are taken over the distribution n(s). If a percolation 
transition occurs at some value p = p c , then in the infinite volume limit, S 
remains finite for p < p c but diverges (it is proportional to the number of lattice 
sites N) for p > p c due to the presence of an infinite cluster. The mean size can 
also be written as the integral of the pair-connectedness function 

1 N 

the quantity within the bra-ket is 1 when i and j lie in the same cluster and 
otherwise. 

One can also consider bond percolation, in which bonds of the lattice are 
randomly colored black with probability p and clusters are then constructed 
from sites connected by black bonds. Bond percolation is generically in the same 
universality class as site percolation. A bond percolation problem on a lattice L, 
for instance, can be mapped exactly to site percolation on the covering lattice 
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of L (as defined in reference ([19])). The properties of bond percolation clusters 
can be deri ved by appropriate averaging with the partition function 

z= £ p b (i-p) Nb -\ N ° (2) 

colorings 

with q — > 1; b equals the number of black bonds that constitute N c clusters in a 
lattice with Nb total bonds. In fact, the partition function of the g-state Potts 
model on this lattice can be recast as (2); thus percolation is a non- interacting 
limit of these theories. The q — > 1 Potts model coupled to gravity has been 
studied analytically as a matrix model by Kazakov [20] . This mapping provides 
an interpretation of typical percolative properties in the language of spin models. 
For instance, the relation (1) is analogous to the equivalence of the susceptibility 
(S) and the integral of the spin-spin correlation function. The scaling behavior 
of this correlation function can then be computed exactly in 2 dimensions via a 
mapping of the Potts model to the Coulomb gas [21]. The result is 

(Sc itCj ) ~ \?i ~ r* r^^+^i), (3) 

with A» >g=1 = A° g=1 = 5/96 \ 

The random surfaces which we shall examine are triangulations of two- 
dimensional geometries. The discretization we use consists of a sum over all pos- 
sible triangulations of N vertices, excluding degenerate triangles 2 , with weights 
determined by the partition function 

Zn = ^ P( T ) ex P(~'Krnatter); (4) 

Umatter is the Hamiltonian for Ising or Gaussian fields and p(T) is the measure 
in the space of triangulations T/v of N vertices. Most of the relevant theoretical 
calculations are performed in the grand-canonical ensemble, with the partition 
function 

oo 

Z(ji) = Z N exp(-fiN) (5) 

JV=1 

dependent on fj,, the cosmological constant. The integrated pair-connectedness 
correlation function in the grand-canonical ensemble then satisifies 

oo 

(X>^>(MW - E NS N Z N cxp(-fiN). (6) 

i,j N=l 
lr The weight A° 1 is not identified with a local operator in percolation theory; it just 
determines the scaling behavior of a certain class of correlation functions. 

2 i.e. those with loops of length 1 or 2 or vertices with fewer than 3 neighbors 
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We will thus simulate the tensor product of the q = 1 Potts model with the 
matter theory coupled to gravity. Only, since perception is not dynamical, it 
does not generate any back-reaction on the gravity orTtie matter. Using standard 
arguments 3 , ([22]) one can express the scaling behavior of the integrated pair- 
connectedness correlation function of percolation coupled to gravity via the KPZ 
formula ([18]). The weight A° x = A° x = 5/96 is dressed by gravity in a 
theory of central charge c = c matter < 1 (c perco ; atJO „ = 0) so that 

\,=i - jjg = U + ^x/T^VI - ^ - ^25=^) J A CT>8=1 (1 - A a , 9=1 ). (7) 

This dressed weight determines the scaling of the integrated pair-connectnedness 
function with fi on surfaces of genus h: 

£> i)C ,)(MW ~ (^-^J^+^^O+^Jd-fc) (8) 

i,3 

(the —1 before the dressed weight accounts for the integrations of i and j over 
the surface) with 

1 



7s= 12 



j-l-V(25-c)(l-c)J. (9) 
Then the relation (6) and 

Z w ~ Ar 1 +(T*- 2 X 1 ->0 (10) 

yield the finite-size scaling relation 

S N ~JV 1 - 2A "«= 1 ; (11) 

this is the scaling law that we shall verify numerically. As more generally in 
spin models, the above exponent is referred to as ^ jvdn. The mean cluster size 
satisfies S ~ (p — p c ) -7 , the correlation length (governed by the decay of the 
pair-connectedness function) £ ~ (p ~ Pc) u and is the intrinsic Hausdorff 
dimension of the random surface being considered. We also shall measure the 
fractal dimension of the largest cluster at p c ; the average maximal size cluster 
of each configuration M. scales as 

M-N^S. (12) 
Standard scaling arguments [24] relate j/i'dn — 2d f J dp — 1, 

' 5 applied to the integrated two point function rather than the one point function as usually 
presented 
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2.1 c = l 



The scaling relations become more complicated for c — 1. The analytic solutions 
of the c = 1 matrix models [23] and a careful analysis of Liouville theory show 
that correlation functions no longer scale simply as powers of the cosmological 
constant /i. Instead, the appropriate scaling variable is r\ which satisfies 4 

[i = — 77 ln(?7) + c\r\ H (13) 

in the limit of small 77; ci is a constant that we do not specify. The analytic 
solution employs a modified Gaussian propagator on phi-cubed lattices and is 
thus not based on the same discretization that we use in our simulations. We 
shall assume that the asymptotic scaling relation also holds for the model we 
simulate; we will not presume universality for the subleading coefficient c\. We 
therefore conjecture that the scaling relation (8) should be modified so that 

£> 4lC ,)(^(M) - ^(-i+A^ l)+(2 - 7s )(i-,). (14) 

Our simulations will be done on worldsheets of toroidal topology (h = 1) for 
which Z(fi) ~ ln(^) for c = 1. To extract the asymptotic scaling behavior 
of Sn, we invert the relation between rj and n order by order in 1/ln/x and 
ln(— ln(fi))/ ln(ju) to obtain 

,. , ln(-ln^) f ln(-\nn) \ 2 ln(-In/x) \ 

! l + + "iw^ + '-j- (15) 

We then expand the inverse Laplace transform of (14) to obtain 

1 / lnlniV In In AT 2 lnlniV 

NS N Z N ~ N{NlnNr [I + ( 1 ^ r ) ~q^2+- 

, u^(-uj) , . In In AT \ . „, 

1+ ^hr-^-^WW + -)^ (16) 

oj = 2(— 1 + A cr . 9= i) = —^/vdn — 1 and 'J is the digamma functiom^The scaling 
behavior of NZn is obtained by inverse Laplace transforming dZ(rj(fi))/d(i: 

1 , 1 lnlnAf . , . 

^^-AT^+hTlV-anW^' - 

In addition to the higher order terms that we have dropped from the inversion, 
there are additional corrections to the above formulae. The corrections to the 




Without loss of generality, we set )i c = 0. 
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logarithmic rcnormalization of /i which depend on c\ in (13) should lead to 
contributions to (16) and (17) that are competitive with the smallest corrections 
we have shown. We neglect the usual corrections to scaling that are suppressed 
by (non-integer) powers of l/N. To compare this finite-size scaling form with 
our data, we shall present a plot of ln(«S2Ar/<Szv)/ln2 computed from (16) and 
(17). The leading logarithmic contributions that we have derived above will turn 
out to be very significant 5 . 



2.2 c>l 

We are particularly interested in using percolation to uncover the character of 
the internal geometry of surfaces for c > 1. In this regard, we discuss the 
behavior that we would anticipate in the large c limit. Ambj0rn, Durhuus, 
Frohlich and Orland [26] have argued that in the large c (d) limit, the string 
path integral (with appropriate measure) should be dominated by surfaces that 
form a branched-polymer like structure. To show this, they essentially perform 
a saddle-point evaluation of the string path integral |~ | 

z= £ (Q*) d/2+tt / LKf ex p -E^EW-^) 2 . as) 

with the integration over the center of mass of X implicitly omitted. T n rep- 
resents the space of triangulations that do not contain loops of unit length or 
vertices of order less than 3. As usual, Cy is the adjacency matrix of the triangu- 
lation T. For i ^ j, it equals the number of links connecting i and j; Ca = — g,; 
qi is the coordination number of site i, equal to the number of links emanat- 
ing from i. We can integrate out the Gaussian fields, obtaining (det' CVj) _d / 2 
(with the zero- mode omitted). It can be shown [26] that this determinant is 
minimized by surfaces consisting of tetrahedra glued together along a common 
pair of vertices. 6 . To perform the gluing, the link connecting these vertices 
within each tetrahedron is slit into two edges; each edge is then glued to a slit- 
ted edge of the adjacent tetrahedron. This process creates loops of length 2 
which are allowed in l~ n , these are excluded ffo3i the set o triangulations T n 
that we consider in our simulations. It has been shown that in d = and —2, 
the critical properties are independent of whether loops of length 2 are included. 
Let us assume, as in reference ([26]) that this [Is~|also true in the large d limit, 
since the configurations that dominate the saddle point in the space T n are quite 

5 The essential role of logarithmic corrections in interpreting numerical measurements of 
7 S at c = 1 has been previously discussed in references [25, 5]. 

6 This does not rule out the possibility of other saddle point solutions, though. 
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Figure 1: A sample triangulated surface built from linked tetrahedra, with 
vertices colored as in percolation. The arcs represent slitted edges. To the right, 
the associated branched polymer appears; each tetrahedron corresponds to a 
polymer link. 



complicated and not amenable to the following simple analysis. We shall then 
analyze a class of triangulations discussed by Ambj0rn et al. [12] that includes 
these saddle point solutions; this class is formed by linearly joining the linked 
tetrahedra that minimize the determinant. A typical triangulation is depicted 
in figure 1. 

We first show that the percolatiEffl-lproblem on this class of trifiHgulations 
can be mapped to a percolation problem on the model of branchc U it olymers 
considered by Ambj0rn et al.. Each tetrahedron will correspond to a bond in 
a branched polymer. We identify an internal polymer vertex with a pair of 
sites connected by a slitted edge on the original triangulation. Vertices on the 
boundary of the polymer are associated with pairs of sites that lie adjacent at 
dead-ends of the triangulation. Now, if either of the two siitesi along a slit edge 
are colored black, then the percolation cluster to which they belong must include 
any of the black sites in adjacent tetrahedra. If an infinite percolation cluster is 
to wind through the lattice, its extent depends only on whether or not at least 
~J one of two sites on the original triangulation, associated with a polymer vertex, 
is black. Therefore, we identify black-black, black-white and white-black pairs 
of sites with a black polymer vertex; a wlrrtg-white pair is mapped to a white 
polymer vertex. If we color sites on the triangulation with probability p* , then 
at least one of a pair of sites is black with probability p(p*) = 1 — (1 — p*) 2 ■ 
Thus if the percolation threshold on the original lattice is p* , then the associated 
branched polymer will undergo a transition at p(p*). 
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We now show that these branched polymers (in their critical limit) do not 
admit a percolation transition for p < 1. They are described by the partition 
function [26] 

2 = EE flX-Nfa)) exp(-MiV), (19) 

n tetjv 

in which fjv represents the set of trees (graphs with no loops) of N vertices 
and each vertex with coordination number qi is given a weight w(qi) = riife)" 
(since the determinant of the adjacency matrix on these surfaces just effectively 
shifts the exponent in the measure). 

The analysis of percolation on these polymers proceeds as in the case of the 
Bethe lattice [24]. First, we randomly select a site i on the lattice and a branch 
that emanates from it. Let T denote the mean number of sites in this branch 
that belong to the same cluster as i. If the neighbor j to site i in this branch 
has coordination number qj then 

T = p + p( gj -l)T; (20) 

the first contribution p is the probability that site i is black and in a cluster; if 
site j is black (with probability p) then each of its qj — 1 subbranches contributes 
on average T sites. Therefore, T = p/(l — (q— If we then randomly select a 
black site on the branched polymer, we see that on average it is part of a cluster 
of size 

g = 1 + <g> r = 1+ (i-(g-i)r) ' (21) 

where S is the mean cluster size. Hence percolation will occur on these trees as 
long as the density of dead ends (q = 1 vertices) is not large enough to reduce 
(q) below 2 and trap all percolating clusters. If (q) > 2, then the mean cluster 
size behaves as if we were just considering a Bethe lattice. 

Now we shall compute the mean coordination number in the critical limit 
for this model of non-interacting branched polymers. First, note that a polymer 
vertex of coordination number q is associated with two sites of the original tri- 
angulated surface with coordination number 3q. If this surface is of infinite size 
and of finite genus, then simply (3q) = 6. One can also show in a slightly more 
painstaking fashion that this condition holds in the grand-canonical ensemble 
by analyzing the solution of the Schwinger-Dyson equations for Z analyzed by 
Ambj0rn et al. [12]. The equations are written down for 'rooted' branch poly- 
mers, which emanate from a marked point with coordination number 1, they are 
represented graphically in figure 2. Z obeys 



exp(-/i)(l + /(Z)); f(Z) =^w(i)Z*-\ (22) 

i=2 
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Figure 2: The Schwinger-Dyson equations for the partition function of a rooted 
branched polymer. 

□ 

Criticality occurs at ii = p, c such that dZ/dii\^ c — > oo providing that / is 
convergent in [0, Z(ll c ) = Z c ]; this criterion is satisfied for the measures that we 
will consider. Z(/i) will then have a leading singularity of the form ~ (/i — /ic) 1 / 2 . 
The mean coordination number (minus 1) equals 

{q } ~ -dz/d» ■ [ 6) 

By differentiating the Schwinger-Dyson equation with respect to w(i) and p and 
substituting in above, we find that simply (q—1) = f'(Z) exp(— li) which equals 
1 — 1/(A). (A) is the mean surface area which diverges at p = p c . Therefore 
S(fi) (21) is non-singular for ll > fi c and p < 1. Thus criticality drives the 
percolation threshold p c to 1. The partition function for the g-state Potts model 
on a branched polymer can also be analyzed in the limit q — ► 1 (this was done 
for the Ising model in ([12])). From this one can obtain an exact formula for the 
number of clusters per unit area, which turns out to just depend linearly on p. 

We close this section by reiterating that the equivalence of the behavior 
of percolation clusters on branched polymers (which we have analyzed) and on 
worldsheets in the large d limit relies on unproven assumptions. It is not perhaps 
so obvious that percolation on the saddle points in the space T n will behave as on 
the saddle points of T ni which we may not even have completely characterized. 
Thus, these large d results should be interpreted with approprfSe caution. The 
preceding arguments do, in any case, illustrate the behavior of percolation on a 
class of triangulations which form tree-like structures. 



3 p c = .5 [— | 

We now turn to an examination of how the geometry of the lattice constrains 
the value of the percolation threshold. Consider the function K(p) = N c /N, the 
number of clusters per unit area of the lattice. Sykes and Essam demonstrated 
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[27] that on any planar triangulated multiply-connected lattice, 

K(p) = K(l-p) + 2p(p- l -){p-l). (24) 

This relation is exact for finite N for planar lattices; for lattices with the topology 
of the torus, it is obeyed up to corrections of order 1/N . The derivation of (24) 
is based on an application of Euler's theorem to the two subgraphs consisting 
of the set of all black (white) clusters. The 1/N corrections just arise from the 
genus term in Euler's equation for subgraphs embedded on surfaces of higher 
topology. The key element of the proof of (24) is the following simple fact: on 
any triangulation, each face consists of vertices which are all connected to each 
other. 

The identity (24) is analogous to the duality relation for the free energy of 
the 2 — d Ising model, though in this case the geometrical notion of lattice duality 
is quite different. We briefly digress to discuss this now. Consider a lattice £ 
and select a set of faces T. We then close-pack all faces T\ cGj3e-packing consists 
of linking all vertices on a face together. This yields a new graph C. Let C* 
then be the lattice constructed by close-packing all faces in C that are not in T . 
Then (24) is a special case of the identity 

K c (p)=K^(l-p) + 4>(p) (25) 

where <f> is referred to as the matching polynomial. C and C* are matching 
lattices; any triangulated lattice is self-matching. 

The analogy with Ising duality has a bit more content, since if there is a 
transition at p = p c , K(p) ~ (p — p c ) 2 ~ a . In the mapping of percolation to the 
q-state Potts model, a is the exponent that describes the singular behavior of 
the free energy (which obeys the Ising duality relation) and the specific heat. 
Thus, we can compute a for percolation coupled to c < 1 gravity with the KPZ 
formula as in (7), substituting the weight of the energy operator (A° = 5/8) in 
the q — ► 1 limit of the Potts model for the spin operator weight of 5/96. Then 
assuming the standard hyperscaling relation, one finds that the dressed weight of 
the energy operator A e equals (1 — a)/(2 — a). For percolation with no matter 
(c = 0), the form of K(p) has been computed by Kazakov [20] by analyzing 
matrix model saddle point equations; a = —2 and K{p) also has logarithmic 
contributions to scaling. From the KPZ formula, we see that the exponent a 
then decreases with c, reaching an irrational value of — (2\/T0 + 2) /3 ~ —2.77 ad 
c= 1. 

Thus for < c < 1, K(p) is always non-analytic at its critical point, and 
perforce by (24), p c = 1/2. If the surfaces admit a percolation transition for 
c > 1 with a value of a that varies continuously with c, then p c must still equal 
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1/2. If we discover that a transition no longer occurs near p c — 1/2, then we can 
conclude that either that there is no percolation transition for < p < 1 (as for 
branched polymers) or possibly that a sort of mean-field percolation transition 
occurs. In the former case, two types of lattice behavior could occur. The first 
possibility would be a collapsed lattice in which virtually all vertices are con- 
nected to each other; such a lattice would have a finite probability of containing, 
an infinite cluster for all p > 0. On the other hand, if the average connectivity 
were sufficiently low (as for the finite-genus surfaces constructed from gluing 
tetrahedra, discussed in the previous section), then no infinite clusters would 
occur for all p < 1 7 . A 'mean-field' transition explicitly reficsJa transition in 
which a is constant and integral over a range of central charge. For instance, in 
flat space, the mean field solution for percolation (above d — 6) yields a = — 1; 
such a solution would describe trees with sufficiently high mean connectivity 
8 . It is not so clear that one can construct a finite genus surface surface in 
the universality class of the Bethe lattice with p < 1; perhaps the mean-field 
scenario is also not allowed for 1/2 < p < 1. In any case, both the mean-field 
and no-transition scenarios do seem to be characteristic of percolation on trees 
rather than surfaces; they should signify a qualitative change in the geometry of 
the worldsheet. 

With this background, we now summarize our plan of attack. We will 
simulate percolation on c < 1 worldsheets, checkijigjthat indeed p c — 1/2 and 
verifying theoretical predictions for critical behavior. Assuming p c — 1/2 will 
allow us to measure the exponents much more precisely than would usually be 
the case. For generically a slight mismeasurement in the critical temperature 
will induce a significant error in the finite-size scaling determination of critical 
exponents. Then to look for a degeneration of the worldsheet, we will attempt 
to determine if p c is no longer near 1/2 for c > 1 . 



4 The Simulation 

We now provide a discussion of details of our simulation and analysis. Our sim- 
ulations were performed on dynamically triangulated random surfaces (DTRS), 
using the standard flip algorithm [29] to update the adjacency matrix Cy. The 
flip algorithm rejected updates that led to triangulations with loops of length 1 
or 2 or vertices with coordination number less than 3. We examined DTRS with 

7 An example of an extremely branched triangulated lattice with p c = 1 also is given in 
reference [28]. | 1 | . 

8 A somewhat pathological example of a triangulated surface in this universalitjd class with 
p c < 1 consists of a lattice of connected tetrahedra with infinite genus (with (q) > 6 and 
Pc = 3/ ( {<?i — 13)) described by the construction in the previous section. The identity (24) is 
no longer (cud for infinite genus surfaces, however. 
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no embedding coordinates as well as with Ising spins or Gaussian fields attached 
to the vertices (with interactions along the links). 

Our action for the Ising model was S — C Yln^ &icrj- We chose a measure 
independent of qi, as in the matrix model formulation. The Ising spins were 
updated using the Swendsen-Wang cluster algorithm [30] . Each update consisted 
of 3N attempted flips of randomly chosen links followed by a Swendsen-Wang 
update of the entire lattice. The critical temperature for Ising spins on lattices 
dual to our triangulations was computed analytically in [31], so by the standard 
Ising duality we can determine that C = (1/2) ln(131/85) ~ .216273. This 
particular discretization of the Ising model has been studied in [32, 17]. 

For c = 1 and 2 we used a standard Metropolis algorithm, in which the 
Gaussian fields were updated via random shifts from a flat distribution. In the 
Gaussian case, 3N randomly chosen flips were performed and then 3N randomly 
chosen X coordinates were updated. We used the standard Gaussian action with 
the conformal measure, which corresponds to the choice a = in (18). We also 
simulated the d = 1 Gaussian model with the choice a = —1/2 and found that 
the critical behavior (as determined from the calculation of finite-size scaling 
exponents) was identical, within our statistics, to that measured for a 0. 

A percolation measurement began with a coloring of the lattice ^wrl h ran- 
dom numbers, produced by a Fibonacci generator. All clusters on lattice 
were then constructed (between adjacent sites with random numbers valued less 
than p) using the method of Hoshen and Kopelman [33] . For smaller lattices and 
smaller values of the central charge, percolation clusters were ©wilt after each 
update of the lattice and matter. For c > 1 we found tr[at|it was more efficient 
to construct clusters after every 5-10 updates of the lattice. For c < 1, clusters 
were formed using several different values of p for each lattice. For smaller lat- 
tices, the values of these observables could be extrapolated to nearby values of p 
using histogramming techniques [34]. Thus, for c < 1, there are strong correla- 
tions between the values of observables at different values of p. Histogramming 
proved to be less reliable for larger lattices and for c > 1. For c > 1, our statis- 
tics were not quite so good and we were thus more wary of statistical artifacts 
influencing our interpretation of the data. We therefore performed entirely in- 
dependent simulations for each value of p, so data points at different p would be 
uncorrelated. 

Our errors were computed using the jacknife technique. We first used the 
jacknife error to estimate the auto-correlation time, r ~ 1/2 {a jack I o naive) 1 \ 
Cnaive is just the standard error, assuming no correlations. Additionally, we 
measured the auto-correlation function, and using standard methods [35] and 
appropriate fits, determined the integrated autocorrelation times (with errors 
on these times obtained via jacknife). Generally, these two techniques yielded 
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consistent results. For some of the largest lattices, the fits to the auto-correlation 
function were more reliable, since the jacknife err or d id not plateau convincingly 
as a function of increasing bin-size. Thus, some l ofl the estimates of the error 
bars for data on the largest lattice (TV = 16384) might be a bit off. 

To obtain the desired precision for our observables, we sampled a large 
number of lattices - between 3 x 10 4 and 2.5 x 10 5 independent samples, requiring 
up to 2 x 10 6 sweeps per data point. We performed the simulations on HP 
9000, IBM RS6000 and DEC MIPS workstations; iir-ttptal we used roughly the 
equivalent of 5 months of time on a HP workstatiorTto gather our data. Our 
lattices were of toroidal topology. For c = 0, we used lattices of size 512 through 
16384; for c = 1/2 our lattices ranged in size from N = 2048 to 16384. The 
correlation times were much longer for c — 1 and 2; thus the c — I runs were 
performed for N = 1024 through 8192 at p c = 1/2 and 1024 and 2048 for other 
p. c = 2 runs were only done on lattices of size 1024 and 2048 and also at 
p= 1/2 using N = 4096. 

Our basic set of observables consisted of : the number of clusters per unit 
volume K(p), the maximum cluster size per configuration A / J(]^)-jind the mean 
cluster size S(p). 

I I Our other main task was to determine the critical temperature, to verify 
that p c = 1/2 for c < 1 and then to look for p c for c > 1. Finding an appropriate 
technique to do this was not so straightforward. To illustrate this difficulty, 
we first recall that the width of the transition region is of magnitude N x / udH . 
This width should be, for instance, roughly comparable to the breadth of the 
susceptibility peab-(%he mean cluster size with the infinite, i.e. largest, cluster 
excluded) and its cTsymptotic shift from p c (N = oo). For c < 1, vdn ranges 
from 4 to about 4.7, indicating that the transition will be extremely broad; for 
lattices of several thousand nodes, dp ~ .1 to .2. In contrast, vd — 2 for the 
2d (flat space) Ising model, which will exhibit a transition that is an order of 
magnitude sharper on lattices of comparable size. 

In fact, significant effort has b een devoted to developing methods to deter- 
mine percolation thresholds precisely] Unfortunately, nearly all of the standard 
techniques are inapplicable in our case. They involve determining the fraction of 
spanning clusters (see e.g. [36]) or comparing the outside and inside perimeters 
of clusters [37] . These methods rely on the identification of a lattice boundary, 
which we cannot do easily on random lattices. Other techniques also exploit 
exact analytic results for small cluster sizes [38], which are not available to us 
in this case. We therefore ended up adopting a variant of a standard technique 
used in spin models, which relies on a determination of intersections of Binder's 
cumulants. 
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and 



Explicitly, wc shall consider the quantities 

= 1 - (26) 

u "^-wm?- (27) 

These quantities should obey (for c < 1 at least) a finite-size scaling form 

usM N >P)~fs,M(z)i z = (p- Vc )N x l vd ». (28) 

/ is the universal finite-size scaling function and z is the appropriate scaling 
variable. The deviation of these cumulants from 2/3 just characterizes the mag- 
nitude of fluctuations of the mean and maximum sizes. For p — > oo, / — > 2/3, 
since in this phase, the infinite cluster saturates the system and fluctuations 
diminish. Asymptotically, the cumulants for different values of N will intersect 
at p = p c (which is much better than p ~ p c ± N 1 / 1 " 1 * 1 ) where they will be linear 
and have slope f^N 1 / 1 " 1 ". Therefore, the slopes will not be particularly large, 
and they will grow slowly; determining the points of intersection will require 
very good statistics. 

We briefly digress to point out how the cumulant us is related to the stan- 
dard Binder's cumulant in spin models, u = 1 — (M 4 )/3((M 2 )) 2 ; M is the 
magnetization. Although it is true that (S) — (M 2 ), in general uj^ us- In the 
Ising model (and this should presumably generalize to all g-state Potts models), 
one can easily derive the relation: 

<M 4 ) = 3<5 2 ) - 1( J n(s) s y J n(s)s); (29) 

n(s) is again the cluster distribution function and the latter term is proportional 
to the probability that 4 black points in the lattice lie together in the same 
cluster. The sum of the terms, which both scale the same way in N, is thus just 
the integrated 4-point function of the spin-operator, with a value that should in 
principle be calculable analytically. So our cumulant just measures one (non- 
negligible) piece of this 4-point function. 



5 Numerical Results 
5.1 Correlation Times 

We now turn to the results of our simulations. First, we remark that we observed 
considerable critical slowing down for A4(p) and S(p) for all values of c. We 
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c = 


c = 1/2 


c = 1 


c = 2 


z M /d H 


.69 (2) 


.75 (6) 


1.5 (2) 


1.4 (2) 


z s /d H 


.70 (2) 


.74 (6) 


1.4 (2) 


1.4 (2) 



Table 1: A summary of the critical slowing down exponent z/dg for our runs 



found that the values of the auto-correlation exponent z/djj '■ t ~ N z / dn 
were identical, within our statistics, for S and M. In the following table we 
summarize the values of these exponents. 

For c = 0, the dynamics consists solely of the local link flip updates. Since 
our cluster observables exhibit considerable slowing down, they reflect large- 
scale correlations in the geometry of the worldsheet. In contrast, measurements 
of the distribution of curvature (e.g. of \qi — 6| ) do not exhibit measurable 
critical slowing down, dn has been estimated to be about 2.87 for c = [40, 3]; 
this leads to a value of z of 2.0. The value z/du ~ -7 is not the largest value 
that one might anticipate. For, one would naively guess that the decorrelation 
of the slowest critical modes of the surface geometry would be determined by 
the spectral dimension, d s — 2 [3], which characterizes the probability that 
a random walker returns to its original site after time t, P(t) ~ t~ da l 2 . Thus 
effectively a random update would take a time proportional to N (corresponding 
to z/du = 1) to diffuse across the surface. 

Critical slowing down becomes more severe with increasing c. This could be 
partially due to changes in the fractal properties of the geometry; presumably 
the updates of the matter induce additional slowing down also. Since the change 
in z/dn from c = 0toc=l/2is rather small, it appears that the Swendsen- 
Wang algorithm is quite efficient for our measurements. Further discussion of 
the efficiency of spin algorithms on random lattices will appear in [41]. Our 
auto-correlation times range from about Ti nt = 1 (for c = 0, N = 512) to 50 (for 
c = 2, N = 2048). Although z/dn is virtually identical for c = 1 and c = 2, the 
correlation times for c = 2 are consistently about a factor of 3 higher than the 
corresponding c — 1 times. Probably, cluster methods would have been more 
efficient in updating the Gaussian fields [42] . 

We observed one other curious difference between auto-correlation functions 
for c = and c > 1. For c = 0, the fits of the auto-correlation function (as a 
function of time T) to the form exp(— T/t) were excellent for T > 3r and all 
values of N. For c > 1, this was no longer true. The auto-correlation functions 
for c = 2 were quite difficult to fit to; even for times considerably larger than 
r, they exhibited very strong transients. We shall argue in subsequent sections 
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that for c > 1 that we do not observe good scaling behavior; the behavior of the 
auto-correlation functions could perhaps be interpreted as a first hint that this 
might be so. 



5.2 Critical Exponents and KPZ 

We now present tables of the effective critical exponents (df /dii)ef f = 
ln(7W2Ar/A^7v)/hi2 followed by (j/vd^eff = hi(6>2JV /Sn)/ In 2 obtained 
through finite-size scaling at p — 1/2, as defined in section 4. 

The exact values for j/udji) obtained from (7) and (11) are 3/4,(4 — 
y/7/2)/3 and 1 — ^/5/24 for c = 0,1/2 and 1 respectively. We have again 
used df/dji = 1/2(1 + ^Ivdn). First, we remark that the agreement between 
the theoretical values and our measurements for c < 1 is surprisingly good, given 
that our runs are on medium-sized lattices. The exponents generally decrease 
as a function of increasing lattice size and central charge c. This seems reason- 
able, since embryonic universes with necks of the size of the cutoff become more 
dominant as c increases; we can see from the theoretical calculations that this 
corresponds to a reduction in the critical exponents. For TV not so much larger 
than the cutoff, the effect of these bottlenecks should be not so apparent and 
then presumably the critical exponents would be over-estimated. For Cp^iO on 
the larger lattices, the exponents have more or less already leveled off,Uimhin 
our statistics, to their asymptotic values. 

For c = 1, if we simply assumed that S and M. scaled as powers of N, then 
our estimates of the critical exponents will be very far off. In fact, the mEaSured 
values of the effective c = 1 exponents agree quite well with the theoretical 
predictions for c = 1/2. The logarithmic corrections, though, will shift the 
effective exponents far away from their asymptotic values. Using the relations 





c = 


c = 1/2 


c = 1 


c = 2 


N = 512 


.879 (2) 




.860 (5) 




N = 1024 


.880 (3) 




.861 (5) 


.81 (1) 


N = 2048 


.880 (3) 


.866 (5) 


.859 (7) 


.83 (2) 


N = 4096 


.871 (4) 


.866 (6) 


.841 (8) 




N = 8192 


.875 (5) 


.852 (8) 






theory 


.875 


.855 


.772 





Table 2: Measurements of (df /du)ef / compared with theoretical predictions 
of the N — > oo limit 
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c = 


c = 1/2 


c = 1 


c = 2 


N = 512 


.774 (2) 




.744 (6) 




N = 1024 


.772 (4) 




.739 (7) 


.68 (2) 


N = 2048 


.768 (5) 


.747 (6) 


.732 (7) 


.67 (3) 


N = 4096 


.753 (5) 


.741 (7) 


.700 (10) 




N = 8192 


.755 (7) 


.722 (10) 






theory 


.75 


.710 


.544 





Table 3: Measurements of ('j/vd^eff compared with theoretical predictions 
of the iV — ► oo limit 



(16) and (17) we see that the first subleading contribution yields 

(<y/vd B ) e „ = 1 /yd H + + • • • . (30) 

In the following plot, we compare our data for (^y/i/d[j)eff at c = 1 with the 
theoretical prediction, derived from (16) and (17). We plot (30) with the dot- 
ted curve. The solid curve takes into account all corrections in (16) and (17); 
the smallest correction in the expansion we used to generate this theoretical 
prediction was of order .002 for N > 1000. At this order further corrections 
arising from the logarithmic renormalization of /i, which we did not compute or 
conjecture to be universal, should appear. We therefore would guess (without 
any rigorous justification) that these additional corrections should probably shift 
the curve by less than the statistical accuracy of our data, which is roughly .01. 
For c < 1 we observed small systematic discrepancies from asymptotic scaling, 
presumably from power law corrections. Since power law corrections have also 
not been taken into account at c = 1, we consider the agreement between our 
measurements and the both c = 1 theoretical curves is thus quite good, certainly 
comparable to the agreement for c < 1 . This suggests the likelihood that at least 
the leading c = 1 logarithmic correction to scaling (30) is indeed universal. We 
emphasize that we cannot predict the presence of logarithms in the finite-size 
scaling ansatz from our data, since it is not nearly precise enough to distinguish 
logarithmic behavior from a small power law. We can clearly see, though, that 
if we assume the logarithmic scaling ansatz given in (16) and (17), that we can 
extract the proper value of j/vdn from our data within an accuracy of a few 
percent. Also, the agreement between our data and the KPZ prediction is quite 
sensitive to the power with which the leading logs appear in the finite-size scaling 
ansatz. 

Naively it appears that since the effective exponents for c = 1/2 and c = 1 
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Figure 3: A comparison of the theoretical prediction of the finite size c = 1 
scaling of with our data. The dotted line includes the leading logarithmic 

correction and the solid line takes into account all subleading terms that we 
calculated. The horizontal dashed line indicates where these curves asymptote. 
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th j.t j re measured are quite similar, that the coupling to matter in both theories 
also is qualitatively of the same magnitude. We see now that this should not 
be tr y^-^ n the continuum limit where the exponents are really very different. 
One would need to simulate^ttle c = 1 model on a 400000 node lattice to see 
{^/vd}{) e fj decrease to .65, which would be a somewhat appreciable, though not 
tremendous difference from the c = 1/2 value of .71; ( n f/i/djj)eff ~ -60 for 10 12 
triangles. One mi glht I think that there is a bit of a conspiracy occurring here: 
the c = 1 KPZ exponents and the logarithmic finite scaling ansatz somehow 
mysteriously combine to produce behavior that strongly resembles c ~ 1/2 on 
lattices that are numerically accessible to us. The discussion in section (6) should 
demystify this observation, though. 



5.3 Determining p c 

We now ascertain if we can successfully estimate p c from our data. First, we 
verify numerically the identity (24), which constrains p c to be 0,1/2 or 1 for 
non-integer a. In figure 4 we have plotted K(l — p) ~ K(p) (recall that K 
equals the number of clusters per site) versus 2p(p— 1/2)(1 — p) using data from 
c = simulations with TV = 8192. Note that all data points in this plot are 
heavily correlated; we used histogramming (from samples taken at p = .49, .50 
and .51) to produce all of the data points. Clearly the agreement is excellent, as 
expected, since (24) should be correct up to corrections of order 1 /N for genus 
1. 

Next, we examine the intersection of the cumulants um{p) and us{p) for 
c = 0, presented in figures 5 and 6 for lattice sizes N = 1024, 2048 and 4096. The 
curves intersect around p c = .506(4) and .508(4) respectively, within a percent 
or two of the exact value of p c — .5. Considering that the large value of vdn 
implies that the width of the transition is of order .15, we conclude that this 
technique is fairly successful in locating the critical point in this case. For 
c = 1/2, this technique is at least moderately successful; the intersection point 
occurs roughly around .51. We demonstrate this in figure 7, which shows the 
curves u M for N = 2048,4096 and 8192. 

At c = 1, the intersection point of us, as shown in figure 8, has moved quite 
a bit to the right. Its position is difficult to determine given our statistics, but 
it appears to be above p = .52 and below .56. The data points in this plot are 
each statistically independent; we have foregone histogramming because it is no 
longer so reliable. Since we know that p c = .5 from equation (24) in this case, 
we can see that the method of intersections now has become much less reliable. 
This failure is expected, though, since the cumulants will only intersect near p c 
when one has reached the scaling regime where (28) is valid. One would expect 
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0.5 0.505 0.51 0.515 0.52 

Probability p 



Figure 4: A comparison of the theoretical prediction relating the the number of 
clusters per unit area to a polynomial in p. The data is from c = simulations 
on lattices of size 8192. 




0.6325 b 



0.49 0.5 0.51 0.52 

Probability p 



Figure 5: The cumulants Um(p) for gravity with no matter. The circles refer to 
lattice size N = 1024, N = 2048 points are labeled by triangles, and the squares 
correspond to N = 4096. 
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Figure 6: The cumulants u s (p) for gravity with no matter. Points are labeled 
as in the preceding figure. 
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0.49 0.5 0.51 0.52 

Probability p 



Figure 7: The cumulants u M (p) for gravity coupled to c = 1/2 matter. Again, 
N — 2048 points are labeled by triangles, TV = 4096 data is represented by 
squares; the circles correspond to N — 8192. 
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0.49 



0.51 0.53 0.55 

Probability p 



0.57 



Figure 8: The cumulants ug(p) for gravity coupled to a Gaussian field. The 
circles and dashed lines represent N — 1024 data and N = 2048 points are 
labeled by triangles connected by dotted lines. Lines are drawn to guide the eye. 
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that the log corrections to scaling would modify (28) to a relation roughly like 



us,M(N,P)~Q- + ^jf)fsMz), (31) 

(the definition of the scaling variable z should also receive logarithmic correc- 
tions). This modification is anticipated because continuum observables (such 
£|S /(0)) seem to generically receive 1/lniV corrections at c = 1. Note that 
as c increases and bottlenecks on the worldsheet become more prominent, the 
values of the cumulants (at p = .5) decrease, indicating that our observables 
experience greater fluctuations. Similarly, increasing N will make the presence 
of necks more apparent; this should lead to a decrease in a cumulant, implying 
that the constant C in (31) is positive. This will then result in a shift in the 
intersections to larger values of p, as observed. 

For c = 2, the intersection moves even further to the right. We illustrate 
this in figure 9; all points in this figure are again statistically independent. 
In fact, we cannot determine whether the cumulants actually intersect or just 
asymptotically merge, since they have not even crossed over each other within 
our statistics at p = .65. Clearly, though, there is no intersection between p = .50 
and p = .56. At first glance, this appears to be exactly what we were looking 
for. Our original goal was to determine if p c was appreciable far from .5. If this 
were true in the continuum, one would infer that either a transition was absent 
for p < 1 or that perhaps the worldsheet admitted a mean-field percolation 
transition; either scenario would be more characteristic of percolation on trees. 
One might hope that the merging of cumulants closer to .6 as observed would 
just indicate that the percolation correlation length had finally grown to well 
beyond the size of our lattices. 

The above scenario, however, is somewhat implausible. Note that the ef- 
fective exponents, e.g. 7 jvdji do not change dramatically at p = .5 from c = 1 
and c = 2. In general, the c = 2 and c = 1 surfaces do not appear qualitatively 
so different, so one might doubt that the percolation correlation length, which 
is infinite at p = .5 and c — 1 has decreased to roughly the lattice size or below 
for c — 2. We saw that the intersection of the cumulants was shifted because 
of finite-size effects at c = 1. Presumably finite-size effects could also play a 
large role at c = 2, so perhaps the large (at least 15% of p\) shift we see in the 
cumulants is just a finite lattice artifact. We shall argue in the following section 
that this is indeed probably true and that in fact any measurements for c some- 
what greater than 1 will primarily be manifestations of lattice artifacts. Let us 
remark that at least we have observed, via these cumulants, that the behavior of 
c < 1 and c = 2 is qualitatively different. The presence of non-analytic behavior 
of K(p) at p — p c along with good scaling behavior on the lattices we consider 
mandates that the cumulants intersect near p = 1/2, as observed for c < 1. This 
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0.55 0.6 
Probability p 



Figure 9: The cumulants us(p) for gravity coupled to two Gaussian fields, 
and lines are labeled as in the preceding figure. 
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clearly is not the case for c = 2. 



j__p Discussion 

To understand what we are seeing at c = 2, we return to a discussion of the 
scaling corrections for c = 1 . A simple physical description of the origin of these 
corrections can be found in Klebanov's review [16]. We begin by considering the 
c = 1 tachyonic operator of momentum p: 

T( P ) ~ J VW\eMwX + (\P\ ~ 2)0); (32) 

<fi is the Liouville field and g is the reference metric (g = gexp(— (f>)). The 
Gaussian field obeys [43] 

(XX) ~ (lniV) 2 . (33) 

Thus, the finite lattice size effectively imposes an infrared cutoff on spacctimc 
momenta of order 1/lniV. The tachyonic energy satisfies 

Therefore, the minimum ground state energy accessible at c = 1 on a lattice of 
size N should be roughly 

for some constant of order unity C . This correction can be quite large; it con- 
tributes to the logarithmic dressing of correlation functions and thus the logarith- 
mic corrections we have measured that mask the degeneration of the worldsheet. 

Furthermore, one would anticipate that for finite N, the Gaussian model 
coupled to gravity should qualitatively resemble a range of c < 1 models. For 
the c < 1 conformal field theories can be expressed as continuum limits of lattice 
models in which the worldsheet is embedded on Dynkin diagrams [44]. c = 
corresponds to an embedding into A2, c = 1/2 to ^3 and c — » 1 is identified 
with the n — ► 00 limit of A n - For fairly large n, the finite worldsheet lattice size 
imposes an effective cutoff on these diagrams and screens out the presence of all 
but a few Dynkin vertices, producing behavior characteristic of lower c. 

The same mechanism should be in force for c > 1. The extrinsic Hausdorff 
dimension, Djj, may no longer be infinite but it presumably is large. If branched 
polymers do dominate in the large c limit, then Dh(c = 1) = 00 and Djjic S> 
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1) = 4 [2, 26] 9 . Numerical work has indicated that for c of about 2 or 3, that 
Dh is roughly 8—10 [45, 46], though these values are probably not reliable 
10 . For c > 1, one would then expect that the minimum ground state energy 
accessible on a lattice of size N would be roughly 

where again C is a constant of order unity, perhaps with a factor of 2tt thrown 
in. For fixed N, one would only expect to observe a gradual increase in the 
density of bottlenecks that pinch off large baby universes as c increases and 
E min effectively decreases. A degeneration of the worldsheet into a polymer- 
like structure would only become clearly apparent when c reaches a value such 
that (36) is approximately zero. Numerically, it has been observed that this 
degeneration is evident roughly when c is about 10-12 for lattices of order N = 
10 3 . If we assume that at this point Dh — 4 [46], then this determines the 
constant C to be about 6 ~ 2tt, an entirely reasonable value. Of course, these 
should be only interpreted as back of the envelope estimates, designed to show 
that this scenario is qualitatively consistent with numerical observations. 

At c = 2, E mm reaches when N ~ c Dh{c = 2 ">12 Dh[c = 2 ^ 2 . Thus, if either C 
is somewhat greater than 1 (e.g. 6 as in the previous paragraph) or if Dh{c = 2) 
is rather large, one would need extraordinarily large lattices to directly see the 
degeneration of the worldsheet. For instance, C = 6 and Dh(c — 2) = 10 imply 
that N ~ 10 13 ; C = 1 and D H (c = 2) = 10 yield N - 2.5 x 10 5 . If C is not 
larger than 1 and Dh(c = 2) is smaller (near 4, for example) then the tachyonic 
degeneration should be observable on the lattices we consider and finite-size 
effects should not be so dominant. These values, however, would then not be 
consistent with previous much larger measurements of Dh(c = 2). Again, we 
do not advocate taking particular numbers too seriously. It is apparent, though, 
that the true continuum behavior at c = 2 might only become evident at scales 
orders of magnitude larger than those that are amenable to simulation. 

One would then expect also that the percolation correlation length at p = .5 
(where a transition would occur for surfaces) should be comparable to the length 
scale at which the surface degenerates into a Ellnched polymer. If we accept 
the possibility that this length scale is very large, then for c = 2, we should 
not be able to detect the presence of this correlation length on surfaces of a few 
thousand nodes. From the relation (36) it follows that we should observe signif- 
icant deviations from scaling, since the shift in the effective ground state mass 



9 This is the Hausdorff dimension for non-interacffilg branched polymers; for low c, if the 
worldsheet is in a branched polymer phase, intera ptiq ns should still be relevant and there 
would be no compelling reason to believe that Djj = d- 

10 We cannot a priori preclude the possibility that Ujj remains infinite for c > 1. 
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is detectable as a function of N (as at c = 1). In fact, l/ln(10 3 ) < l/(10 3/jD ) 
for D > 4, so for the range of N that we simulate, the shifts in finite-size cor- 
rections could be larger for c = 2 than for c = 1. Therefore, the large shift in 
the cumulant intersections at c = 2 is most likely just a manifestation of a lack 
of good scaling behavior. All other possibilities do not fit very well with at least 
some portion of the prevailing numerical evidence. If at c = 2 the worldsheet 
had degenerated into a tree-like structure, and our measurements were indeed 
satisfactorily measuring asymptotic scaling behavior, then some evidence of the 
tachyonic degeneration should have been easily apparent in previous simulations. 
If on the other hand, no such surface degeneration had occurred and our mea- 
surements were within the scaling regime, then we would have not anticipated 
a large shift in the cumulant intersection from p = .5. 

This leads us to a rather disappointing conclusion. The above arguments 
suggest the likelihood that generically numerical simulations for c somewhat 
greater than one will just measure lattice artifacts. This has also been borne 
out by previous attempts to measure j s , which indicated that corrections to 
scaling were predominant for c > 1 [25]. On lattices of accessible size, the onset 
of the tachyon could be hidden by finite-size effects, and we would fail even to 
capture the qualitative nature of these models in the continuum limit. Scaling 
exponents and critical properties of the geometry (the Hausdorff dimension, e.g.) 
measured in these simulations would not describe continuum scaling behavior. 
For very large c, simulations may reflect the character of the tachyonic instability, 
however. At c = 1, the results of simulations can only be properly understood if 
one knows the form of the corrections to scaling, which can be computed because 
th e theory is solvable. For c < 1, we have found that numerical simulations do 
I — reproduce, with reasonable precision, critical exponents that characterize the 
behavior of percolation on strings. 
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